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Abstract 

Strong nonlinear effects combined with diffusive coupling may give rise to unpredictable evolution in 
spatially extended deterministic dynamical systems even in the presence of a fully negative spectrum 
of Lyapunov exponents. This regime, denoted as "stable chaos" , has been so far mainly characterized 
by numerical studies. In this manuscript we investigate the mechanisms that are at the basis of this 
form of unpredictable evolution generated by a nonlinear information flow through the boundaries. 
In order to clarify how linear stability can coexist with nonlinear instability, we construct a suitable 
stochastic model. In the absence of spatial coupling, the model does not reveal the existence of any 
self-sustained chaotic phase. Nevertheless, already this simple regime reveals peculiar differences 
between the behaviour of finite-size and that of infinitesimal perturbations. A mean-field analysis 
of the truly spatially extended case clarifies that the onset of chaotic behaviour can be traced back 
to the diffusion process that tends to shift the growth rate of finite perturbations from the quenched 
to the annealed average. The possible characterization of the transition as the onset of directed 
percolation is also briefly discussed as well as the connections with a synchronization transition. 
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1 Introduction 



Unpredictable evolution in dynamical systems is due to the propagation of in- 
formation. For instance, the sensitivity of trajectories to infinitesimal pertur- 
bations implies that any arbitrarily small inaccuracy in the determination of 
the initial conditions is exponentially amplified in time, with an average rate 
associated to the positive component of the Lyapunov spectrum. The integral 
of this component, the so-called Kolmogorov-Sinai entropy Q, measures the 
production rate of information that flows from the less significant digits of the 
dynamical variables to the more significant ones. In particular, the existence 
of at least one positive Lyapunov exponent is a sufficient condition for identi- 
fying chaotic dynamics; conversely, a fully negative spectrum is suggestive of a 
periodic evolution. 

This approach to unpredicatable evolution based on linear stability analysis 
has been developed in the context of finite (low-dimensional) systems. Its exten- 
sion to spatially extended dynamical systems is based on the implcit assumption 
that they can be viewed as a collection of almost independent, finite dimensional 
subsystems. The existence of a limit Lyapunov spectrum ^ provides strong 
support to this hypothesis and typical chaoticity indicators, like entropies and 
generalized dimensions, can be turned into their corresponding densities Q]. In 
fact, a primary interest in the study of space-time chaos is the identification of 
thermodynamic- like properties. 

However, Lyapunov instability is not the only source of unpredictability in 
such systems. Actually, information can also flow through the boundaries and 
be transmitted in space by nonlinear mechanisms of front propagation. The so- 
called chaotic rules of Deterministic Cellular Automata (DCA) Q are typical 
examples of unpredictable evolution in the absence of linear instabilities. There, 
the discreteness of the state variable prevents the very existence of infinitesimale 
perturbation, while, on the other hand, even isolated "state-flips" may propa- 
gate through the lattice with finite velocity, giving rise to an irregular dynamics. 
In fact, any DCA rule defined over a lattice of L cells is bounded to exhibit a 
periodic behaviour, since the number of possible states is finite [h^ if h is the 
number of possible states in each given site). What makes a "chaotic" DCA 
rule different from an ordered one is the exponential growth of the recurrence 
time of the typical configurations. Accordingly, the unpredictable behaviour is 
dynamically persistent only in the infinite-size limit. 

A very similar unpredictable behaviour, denoted as stable chaos, has been 
observed also in coupled map lattice (CML) models j^, |7| of the type 

a^^(^ + l) = /(|a:^_l(^) + (I-e)a:,(^) + |a;,+l(^)) , (1) 

where Xi{t) is the continuous state variable at time t on the site « of a Id lat- 
tice; the function / is a linearly stable map of the interval [0, 1] into itself 
and £ is the strength of the diffusive coupling. It is worth pointing out that 
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this spatial coupling cannot produce any linear instability mechanism and the 
whole spectrum of Lyapunov exponents is found to be negative. Accordingly, 
the CML dynamics must eventually approach a periodic stable attractor. This 
notwithstanding, if / is equipped with a sufficiently strong nonlinearity (e.g. a 
discontinuity or a region with rapidly varying slope), one can find a region in 
parameter space, where the "transient" evolution towards the periodic attractor 
grows exponentially with the system size L, analogously to chaotic DCA rules. 
Despite there is no rigorous proof of this statement, many independent numer- 
ical studies confirm such a scenario (see, e.g., ^,0). Striking features of this 
"transient" regime are its stationarity and apparent ergodicity: for instance, 
space and time correlation functions decay exponentially like in usual chaotic 
phases, while ensemble avarages coincide with time averages (provided that 
a size- independent pre-transient is discarded). Moreover, the maximum Lya- 
punov exponent is found to approach a stationary negative value, before quite 
suddenly turning to the value corresponding to the eventual attractor. Because 
of its exponential growth, the "transient" represent the truly relevant regime 
in the thermodynamic limit, while the periodic attractor(s) have no practical 
significance. 

A careful inspection indicates that the mechanism of information production 
in stable chaos is a flow from the outer (left and right) parts of the chain, like in 
chaotic DCA rules [Q. Accordingly, the unpredictability of stable chaos relies 
on a genuine nonlinear propagation mechanism. 

As we recall in section II, stable chaotic evolution can be detected by mea- 
suring damage spreading, i.e. the average velocity of a front propagating into 
an unperturbed region. This indicator can be viewed as a sort of generalization 
of the standard Lyapunov exponent and the use of a proper metric, attributing 
an increasingly smaller weight to the farther sites would make the correspon- 
dence more transparent Q. However, it is worth noting that the front velocity 
does not allow defining an analogous of a negative Lyapunov exponent, since 
in the ordered phase a perturbation does not only regress but also decreases 
everywhere in size leading to the disappearance of the front itself. 

Numerical studies of stable chaos have contributed to shed some light on 
the relationship between this dynamical regime and the appearance of many 
interesting complex phenomena, such as nonequilibrium phase transitions, spi- 
ral chaos, and the propagation of rough interfaces |l^, |ll[. However, little is 
rigorously known about the underlying mechanisms. The aim of this paper is 
precisely to make some progress in this direction by investigating the conditions 
under which finite perturbations can propagate instead of die out. In analogy 
to various analytical techniques introduced to estimate the maximum Lyapunov 
exponent in standard chaotic CMLs |lj, |l^ , here we assume that the evolution 
in phase-space generates a truly random pattern, characterized by short range 
correlations and thereby introduce a suitable stochastic model describing the 
evolution in the difference space. Besides diffusion, the model dynamics allows 
for a random alternancy of a contraction and an expansion process the proba- 
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bility of which depends on the perturbation size. First of all, in Sec, II we verify 
that, with an appropriate choice of the parameter values, the model is able to 
reproduce not only qualitatively but also quantitatively the main features of 
stable chaos. Afterwards, we further simplify the rule determining the expan- 
sion process to reduce the subsequent technicalities, while keeping untouched 
the key ingredients of the model. 

In order to understand how a finite perturbation can be sustained even in the 
presence of an average contraction rate, in Sec. Ill we first discuss the uncoupled, 
i.e. 0-dimensional, case, where, by definition, propagation of perturbations is 
absent. The negative Lyapunov exponent obviously implies that perturbations 
are eventually absorbed, so that stable chaos cannot exist in this framework. 
Nonetheless, the presence of a non-uniform contraction process yields nontriv- 
ial properties of the dynamics. They are exemplified by the difference existing 
between the standard multifractal spectrum (associated to infinitesimal pertur- 
bations) and the spectrum defined in this paper to describe the evolution of 
finite perturbations. In fact, a noteworthy result of Sec. Ill is that it is pos- 
sible to define a finite-size multifractal spectrum independently of the initial 
amplitude of the perturbation. Finally, in Section III we comment about the 
connection with the finite-size Lyapunov exponents recently proposed by some 
authors as a tool to characterize dynamical unpredictability, beyond standard 



linear stability analysis Il4, 15 



The methods and concepts introduced in Section III are applied in Section 
IV to the study of the spatially extended version of the stochastic model. For 
the sake of space, we limit our discussion to the case of "democratic" coupling, 
but it is clear that stable chaos arises in a broad region of strong coupling. 
It is precisely the diffusive coupling to be responsible for the sustainment of 
stable chaos. Over small scales, contraction is more effective than the sporadic 
amplification: in such circumstances, diffusion just levels the damping process. 
Conversely, at larger scales, diffusion proves to be an effective mechamism to 
propagate locally generated amplifications. The net effect is that, in suitable 
parameter regions, perturbations self-sustain. 

By introducing a factorization hypothesis of the spatial degrees of freedom, 
we obtain a good estimate of the probability distribution associated with the 
spatially extended dynamics. In particular, we find a critical value of the con- 
traction rate separating a dynamical regime, where any perturbation eventually 
vanishes, from a truly chaotic phase, where the average value of the perturba- 
tion remains finite in the infinite-time limit. A mean field argument provides 
a suggestive description of this scenario: the diffusive coupling induces a shift 
of the finite-size Lyapunov exponent that grows from the negative quenched 
average (corresponding to the uncoupled limit) towards the annealed average 
of the expansion rate. Depending on the parameter values, this latter average 
may be strictly larger than 0, thus implying that the zero- amplitude regime is 
unstable. 

The transition exibited by the stochastic model is reminiscent of the "fuzzy" 
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Figure 1: Picture of the local deterministic map (g) for typical parameter 
values. 



transition region found in a deterministic CML However, we cannot push 
the analogy to a quantitative level, since the stochastic CML dynamics is self- 
generated and thus it is increasingly regular while approaching the transition 
region (this is the main reason for the fuzziness observed in Ref. ]l7t), while 
here the stochastic properties of the contraction/expansion process is fixed a 
priori. A tight analogy, instead, exists with the stochastic synchronization in- 
duced by an additive noise Indeed, in this context, the external noise does 
not change its robust stochastic features when passing from the synchronized 
to non-synchronized regime. It is precisely this analogy which suggests that 
the transition described in this paper should belong to the universality class of 
directed percolation The numerical simulations described in Sec. V do con- 
firm such an expectation, but subtle problems still prevent us from still drawing 
definite conclusions. Such open problems and the possible future perspectives 
are summarized in Sec. V. 



2 Generalities 



The typical functions used to investigate stable chaos in CML are piecewise 
linear maps of the unit interval of the form (see also Fig. (IT]) for a pictorial 
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representation) 

{cix < x < 1/ci 

l-(a;-l/ci)(l-C2)/A 1/ci < x < 1/ci + A . (2) 

C2 + C3[a;-(l/ci + A)] 1/ci + A < a; < 1 

The most studied case in the hterature corresponds to the hmit A ^ 0, where 
f{x) reduces to a discontinuous function. As we shaU comment along this paper, 
a strong nonhnear component of map f{x), rather than a true discontinuity, is 
sufficient to yield stable chaos. This is why we prefer to consider here the more 
general case . The parameter range of interest for the present study is when 
all initial conditions (except for a set of zero Lebesgue measure) converge to the 
same periodic orbit. For instance, for ci = 2.7, C2 — 0.07, C3 = 0.1 and A = 0.01, 
a stable period-3 orbit exists with Lyapunov exponent A = —0.316... As already 
mentioned in the Introduction, there is a range of values of the diffusive coupling 
£ defined in (|]) , where the CML dynamics exhibits a "chaotic" evolution, despite 
the largest Lyapunov exponent is negative. Damage spreading analysis provides 
a first hint about this mechanism, that is responsible for the sustainment of 
irregular behaviour. More precisely, while standard chaos amounts to a flow of 
information from the less to the more significant digits, stable chaos is generated 
by a flow of information from the outer (left and right) to the inner parts of 
an infinite chain. Unfortunately, while the former flow can be "easily" studied 
thanks to the linearity of the process (in fact, over sufficiently small scales, any 
smooth function can be linearized), the same argument does not apply to the 
dynamics on the left and right edges, that is equally nonlinear at any spatial 
position. 

This is the main reason for the difficulty in deriving the necessary and suf- 
ficient conditions for the propagation of perturbations. Moreover, while the 
dynamics of an infinitesimal perturbation can be studied by neglecting propa- 
gation phenomena, the opposite is not possible. One cannot study propagation 
without properly accounting for the local contraction/expansion mechanisms. 

The damage spreading analysis is performed by studying the dynamical vari- 
able Ui{t) = \xi{t) — yi{t)\, i.e. the absolute value of the difference between two 
test trajectories, Xi{t) and yi(i), that are initially set equal to one another on 
the right of some lattice site, say i = 0, while they are assumed to be totally 
independent on the left. With such a setting, the damage spreading analysis 
amounts to studying the propagation of a front separating the region in space 
where Ui{t) ^ 0{1) (the tail) from the region where the two trajectories converge 
to each other, Ui{t) ^ (the forefront). 

The first conceptual problem that we have to face is not just the propagation 
of the front, but its self-sustainment in spite of the local average contraction rate. 
In order to shed some light on this crucial point, we have simplified the model by 
assuming that the dynamics is indeed irregular, thereby determining whether 
this assumption is consistent with the sustainment of an 0(1) perturbation. 
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This is analogous to the consistency approach developed for the description of 
standard chaos, where complete randomness of the multipliers of the evolution 
operator in the tangent space is assumed in order to estimate the maximum 
Lyapunov exponent |l3| . 

Accordingly, we introduce a suitable stochastic model to describe the evo- 
lution in the "difference" space spanned by Ui(t). Let us start from the simple 
case A = 0. With reference to the CML dynamics, the dynamical rule is com- 
posed of two steps. The first one corresponds to the application of the standard 
discrete diffusive operator 

= (l-£)u,(t) + |iii+i(i) + |u,_i(<) , (3) 

where < e < 1 is the coupling parameter. The second step contains the 
stochastic component of the evolution rule, 

r r.jS , w.p. p^min [1,6^.(0] 

[ aui{t), w.p. 1-p ^ ' 

where, "w.p." is the shorthand notation for "with probability" , while ri{ui, t) is 
a random number distributed in the unit interval according to some probability 
distribution that depends on Ui. 

This is the non-trivial part of the stochastic model, defined to mimick the 
evolution of perturbation in the CML dynamics. The first line describes the 
instability mechanism associated with the discontinuity while the second line 
describes the contraction of ui (t) by a constant factor |^ < a < 1 . In the CML 
model, the instability mechanism arises whenever the test trajectories, Xi(t) 
and yi{t), lie on different sides of the discontinuity of the map. In this case, the 
value taken by the difference variable Ui{t + 1) is not uniquely determined by 
the value of Ui{t), since it depends also on Xi{t). In particular, for small values 
of Ui{t), the instability mechanism occurs quite rarely (it is unlike that the 
discontinuity is placed across two nearby trajectories) and Wi(t-f 1) is amplified 
by a big factor; conversely, for large values of Ui{t), the discontinuity plays a 
role much more frequently, but it is less effective. The numerical analysis of the 
CML model shows that the probability p of the instability mechanism grows 
linearly for small values of the difference variable Ui{t) and approaches 1 when 
Ui{t) does the same. Accordingly, in the stochastic model (^), we have decided 
to schematize this dependence with the simple law p = min[l, bui{t)\. 

Such a stochastic model can be straightforwardly generalized to cover the 
case corresponding to a CML dynamics with a nonzero A. The main difference 
is that, whenever u is smaller than A, u is expanded by a fixed factor A"-'^ 
with a constant probability 6A, while the previous stochastic rule still applies 
for M > A. This amounts to assuming that only perturbations larger than A 



^For the sake of simplicity we assume that the contraction rate is a constant, as in 
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perceive the steep branch of the map as an effective discontinuity. Accordingly, 
Eq. (0) is replaced by 




ri{ui,t), w.p. p = min [I, bui{t)] 

aui{t), w.p. 1 — p 

Ui{t)/A, w.p. p — bA 

aui{t), w.p. 1 — p 



if > A, 

(5) 

if Ui{t) < A. 



One can easily check that these formulae reduce to Eq. (Q) in the limit A — > 0. 

The maximum Lyapunov exponent of the CML dynamics corresponds to 
the time average of the expansion rates of infinitesimal perturbations. In the 
context of the stochastic model, it naturally corresponds to the quantity 

Ao=lim lim lyinMLllM^ lin, lim ^ In M% , (6) 

t~.oo „(0)^0 i ^ \\u,{t)\\ t-oo„(0)^0t ||Uz(0)|| 

that we still denote as the "maximum Lyapunov exponent" . 

An analytical estimate of Aq can be obtained by a mean-field argument, 
according to which the probability of applying the expansion factor 1/A is bA, 
while the probability of appling the contraction factor a is 1 — bA. One easily 
obtains 

Ao«6Aln^^^ +(l-&A)lna = lna-5Aln(aA) . (7) 

Numerical simulations indicate that the true value of Ao is generally slightly 
larger than this mean-field estimate. 

Moreover, for A = 0, Aq « In a < 0, while for increasing values of A it 
may become positive, indicating that a standard choatic regime is attained^ 
Here, we are interested in studying only the parameter region where Aq remains 
negative. 



2.1 Comparison between stochastic and deterministic mod- 
els 

The reliability of the stochastic models (0) and (j|) has been tested by numerical 
simulations which show that both rules (|4|) and (||) together with the diffusive 
coupling (|^), exhibit the same qualitative features of the CML dynamics (^, (H). 
In particular, when the contraction is relatively strong, any initial perturbation 
Ui{0) is quickly absorbed to the fixed point Ui — V«: this regime corresponds 
to the "ordered" or, equivalently, to the "synchronized" phase. On the other 
hand, for weaker contractions, almost any initial condition evolves towards an 

■^This is not entirely correct, since, as discussed in Ref. [po|], the standard chaotic regime 
occurs when a stronger constraint is met: the linear propagation velocity coincides with the 
nonlinear one. 
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Figure 2: Snapshot of the perturbation profile in the chaotic phase of the 
stochastic model ||) at time t — 500. 

irregular spatial structure like the one shown in Fig. ||. In this case, the space 
average 

1 ^ 

u{t)^-Y,Mt) (8) 

remains finite and u{t) = u* , independently of t and L for sufhciently large sizes. 
To be more precise, perturbations of the order 0(1) survive only for a finite time 
r also in the chaotic phase of the stochastic model. Nevertheless, in analogy to 
the CML model, r grows exponentially with the lattice size L. Furthermore, 
ensemble averages indicate that this regime is asymptotically stationary in time 
and the comparison with time averages (obviously performed over times much 
shorter than r) indicate that ergodicity holds as well. It is therefore meaningful 
to define the single-site probability distribution Q{t, v) of finding a perturbation 
Ui in between v and v + dv at time t, and its stationary limit Q{v), attained 
for large enough times. It is also worth introducing the space and ensemble 
averages m{t) = {u{t)) (from here on, unless otherwise stated, (•) denotes an 
ensemble average) and the stationary limit, m(t) — s- to. 

In the numerical investigations we have assumed no-flux boundary conditions 
as they preserve the u = Q fixed point and are reasonably harmless in the chaotic 
phase. The u-dependent probability distribution of the random variable ri{ui, t) 
has been reconstructed from the CML dynamics. Moreover, in order to avoid 
boundary effects in the damage spreading analysis, the lattice size has always 
been chosen in such a way that the perturbation front never reaches the lattice 
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Figure 3: Average front profiles (upper curves) and the corresponding variances 
(lower curves) at three different times, t — 500, t — 1000 and t — 1500. Panel 
(a) refers to the deterministic CML model (P,(^), with parameter values ci = 
2.7, C2 — 0.07, C3 — 0.1, A = 0, e = 2/3; averages have been performed 
over 10^ initial conditions. Panel (b) refers to the stochastic model (^, (^, 
with parameter values a = 0.9, b = 1.7, A = 0, e = 2/3; averages have been 
performed over 10^ realizations, while the probability distribution for ri{ui,t) 
and the values of parameters a and b have been determined from the CML 
model depicted in the left panel. 

edges during each simulation. 

In order to investigate damage spreading phenomena, the initial conditions 
have been fixed by imposing Ui(0) — for i > and randomly choosing Ui{0) 
with a uniform probability distribution in the interval [0, 1] for j < 0. The 
statistical fluctuations showed by Ui{t) have been smoothed out by performing 
ensemble averages (over different realizations of the stochastic process) of the 
spatial configurations at equal times. In the chaotic phase, the initial "kink"- 
like structure persists: the front connecting the perturbed with the synchronized 
region moves with a fluctuating velocity. The ensemble averages corresponding 
to three different times are reported in Fig. ^d, where they are compared with 
the results directly obtained from the CML model (for the sake of space we limit 
ourselves to consider the case A = Oa), where ensemble averages are performed 
over different initial conditions. Several observations are in order. First, we 
notice that, in spite of the simplifications introduced in the stochastic model 
(besides the lack of space-time correlations, we have indeed assumed a constant 
contraction rate as if the two branches of the local map (|) had the same slope). 
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Figure 4: The average front size w{t) of the stochastic (full line) and CML 
(dashed line) models described in the caption of Fig. ^ plotted versus the square 
root of time. See the text for the definition of w{t). 

there is a reasonable agreement with the CML data. In particular, it can be 
seen that the average height u* is definitely smaller than 1 in both cases: the 
reason can be traced back to the combined effect of the contraction mechanism 
with the diffusive process. Somehow larger differences can be observed in the 
behaviour of the variance cr^ — {{ui{t))'^) — {ui{t))'^ but they can be attributed 
more to the approximation in the description of the perturbed region, rather 
than to peculiarities of the propagation. 

In fact, if we look at the width w{t) of the average front, defined as the 
distance between the rightmost sites where {ui{t)) is larger than 0.8u* and, 
respectively, 0.2u*, we observe a nearly square-root growth in both the CML 
and the stochastic model (see Fig. ^, where a nearly quantitative agreement 
is also observed). Moreover, we have verified that, keeping the contraction 
parameter a fixed, and increasing A up to small values (e.g. A ~ 0.01), vp 
slightly increases in both models. 

On the basis of the results reported in Fig. U, one can summarize our obser- 
vations of the front dynamics by effectively assuming a stepwise shape for the 
profile at any time and approximating its motion with a diffusive process with 
drift. More precisely, we can write 

u* /■+°° / \ 
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where vp is the average front velocity, while D is the "diffusion" coefficient 

accounting for the square-root growth in time for the standard deviation of the 
front position x (assumed to be continuous, for the sake of simplicity). 



2.2 Further simplifications of the stochastic models 

As we have verified that the qualitative features of the front dynamics do not 
depend on the shape of the probability distribution of ri{ui, t), we have decided 
to simplify the stochastic models by assuming a S-like distribution, i.e. that 
ri{ui,t) = 1. Moreover, in order to got rid of unnecessary technical complica- 
tions we shall assume that b = a. Notice also that, having chosen niin^t) = 1, 
we are obliged to assume & < 1, in order to avoid the appearence of a fictitious 
fixed point u* = 1. Prom now on we study the following model 



that we call Continuous Stochastic Model (CSM). In the limit A ^ it will be 
named Discontinuous Stochastic Model (DSM). 

The only parameter that we are going to consider in the following sec- 
tions is the contraction parameter a: this is sufficient to identify and char- 
acterize the relevant transition from the chaotic/unsynchronized phase to the 
ordered/synchronized one. 

3 The zero dimensional case 

The most important problem of linearly stable chaos is to understand how fi- 
nite perturbations can propagate in spite of the average local contraction. To 
clarify this point, in this section we consider the zero-coupling limit, i.e. the 0- 
dimensional case. We shall see that even if "chaotic" motion cannot be sustained 
whenever the Lyapunov exponent is negative, the lack of a uniform contraction 
induces anyhow non- trivial properties. 

Our first observation concerns the probability P{t, w(0)) for a finite-amplitude 
perturbation u{0) to be never amplified by the instability mechanism over a time 
t. In the DSM, such a probability can be easily factorized as 



= (1 - s)Ui{t) + -Ui+i{t) + -Wi_i 



it) 




(10) 



P{t, u(0)) = n (1 - "(O)"") = exp ^"(1 " "(O)"") 
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exp 



k=l n=l 

t — »C30 

— > exp 



exp 



E 



k=l ^ 



fe=l 

k 



k 1 — a'^ 



(11) 



One can see that, for t oo, P{t, u{0)) approaches a finite value P(u(Q)) that is 
both strictly larger than and smaller than 1 for any value of u(0) (if < a < 1). 
The same conclusion can be drawn also for the CSM, although the algebra is 
more complicate in that case. The inequality P{u{0)) > indicates that the 
occasional amplifications are not so strong as to prevent the eventual absorption 
of the perturbation (this is consistent with our goal to deal with linearly stable 
processes). On the other hand, the inequality P(m(0)) < 1 indicates that the 
amplification process cannot be neglected. Notice that, since this holds true 
independently of u{0) (although P{u{0)) 1 for u{0) — > 0), there is a differ- 
ence with the convergence to a stable fixed point in a topologically chaotic map 
(think, e.g., of the logistic map in one of the stability windows that follow the 
first period doubling) since, in that case, there would be a threshold (corre- 
sponding to the border of the basin of attraction), below which a monotonous 
contraction would start. A closer similarity exists with the so-called strange 
nonchaotic attractors, as they are characterized by a nonmonotonous contrac- 
tion even arbitrarily close to the attractor ^ . 

The multifractal theory provides the most appropriate framework to 
characterize this system. By following this approach, devised with reference to 
infinitesimal perturbations, we introduce the exponential growth rate of a finite 
initial perturbation u{0), 



(12) 



In the limit u(0) 0, A(t, u(0)) ^{Lh ^^'^ equivalent in this context of the 
usual finite-time Lyapunov exponent ||24| . The proper indicator to look at is 
the probability distribution P(A, t,M(0)) to find a growth rate between A and 
A-|-dA at time t starting from u(0). More precisely, we introduce the finite-size 
multifractal spectrum 



H{A) = lim 



ln7'(A,<,u(0)) 



(13) 



As it will become clear later, H{A) is independent of u{0) (as long as u{0) > 
0) . Nevertheless, we shall show that it differs from the standard multifractal 
spectrum h{X), obtained by taking the limit m(0) before the infinite-time 
limit 

'1 



h{\) = hm 



lim 

oo u(0)^0 



-ln7'(A,t,u(0)) 



(14) 
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In other words, the order of the two hmits is crucial for understanding the 
difference between the behaviour of infinitesimal and finite perturbations. In the 
following, we shall compare the multifractal distribution of X(t) [linear analysis) 
with that of A(t, u(0)) {nonlinear analysis) in both stochastic models. 



3.1 Linear analysis 



We start from the CSM, as the DSM is nothing but a limit case of the former 
one. The linear approximation amounts to assuming Ui{t) < A Vi,<. In this 
case, standard combinatorial analysis implies 



hAiX) 



ln(A) + A 
In(aA) 
ln(a 



ln(l - aA) 



/ln(A) + A 
V In(aA) 



+ 



In(aA) 



In 



aAln(aA) 



ln(a) - A 



(15) 



where we have made explicit the dependence on the parameter A. In the limit 
A ^ 0, the above expression reduces to 



ho{\) = lim - \nV[X{t)] ln(a) - A, A > Ina 

t^oo t 



(16) 



On the other hand, performing the limit A — > before the t oo limit would 
yield the trivial result ho{X) — with the support of ho restricted to the point 
A = Ina. Thus, the non-commutativity of the limits A — > and i — s- oo reveals 
that the discontinuous case is a singular limit of the CSM. In other words, 
the multifractal spectrum of the discontinuous model depends on the way it is 
defined. We prefer to adopt a "physical" point of view, i.e. to consider the 
discontinuity as the limit of a negligible A, which corresponds to taking first 
the i — !■ oo limit. 



3.2 Nonlinear analysis 

We now consider a finite perturbation u{0) in the simple context A = 0. From 
Eq. ( p^ ) we see that A is uniquely determined from the knowledge of u{t) and 
u{0) (and, obviously, the time t). Accordingly, the knowledge of V{A,t,u(0)) 
is fully equivalent to that of the single-site probability distribution Q{t, u), to- 
gether with the initial condition 

Q{0,u)^d{u-u{0)) . (17) 

It proves useful to introduce the notation u = a" and u{0) — a"'% where n and 



no are real variables > 0. Eq. (12) can be rewritten as 



X^X{t,n,no) = ]\n(^^=^^\na . (18) 
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Once Q{t, n) is known, Eq. ( Jlq ) allows one reconstructing the corresponding 
probability distribution V{A,t,u{0)). Two possibilities are in order, either the 
system has never been "kicked" , i.e. reset to 1 by the instability mechanism, in 
which case the initial value has been contracted t-times by a factor a, or it has 
received at least one kick, loosing memory of the initial condition. In the former 
case, occurring with probability P{t,u{0)), A(i, u(0)) = Ina or, equivalently, 
n = no + t. In the latter case, A(t, u(0)) > Ina and the accessible values of n 
are restricted to positive integer numbers strictly smaller than riQ + t (since that 
the maximum possible contraction factor in t time steps is a*). Accordingly, n 
can be interpreted as the elapsed time since the last kick, and the probability 
distribution Q{t,u) can be factorized as the product of the probability G{t — 
n, no) of receiving a kick at time t — n for an initial perturbation u(0) = a"" by 
the probability P(n, 1) of not being kicked anymore for the remaining n time 
steps, 

Q{t,n)^G{t~n,no)P{n,l) (19) 

(in the remaining part of this section and in App. A, with no ambiguity arising, 
we denote with Q also the probability density of the logarithmic variable n). 
Moreover, we impose the condition P(0,u(0)) :~ 1 to extend the validity of 
Eq. (|l9|) to the case n — 0. 

The probability G(t,no) can be recursively expressed as the probability of 
receiving the very first kick at time t plus the probability of receiving the second 
last kick at any previous time, i.e. 

t-i t-i 
G(t,no) =a*+"° Y[{1 - a''+"'>) +J2G{k,no)a'-''P{t - 1 - k,l), (20) 
fc=i fc=i 

with G(l,no) := a^°^^. Eq. (pO[) can be numerically iterated in the large time 



limit to obtain the expression of the multifractal distribution (13). The analysis 
performed in App. ^ shows that H{A), defined as in (|l3|), is a segment of straight 
line restricted to the open interval (Ina. 0) of negative values. The slope of this 
straight line depends on the contraction parameter a but is independent of u(0). 
In fact, in the appendix we show that the slope can be obtained by solving an 
eigenvalue problem, where u{0) enters to specify the initial condition but not 
the operator itself. An approximate analytic solution is also determined, which 
confirms the numerical observation that the slope increases monotonously from 
-1 (for a — > 0) to (for a ^ 1). Curves 2 and 4 in Fig. || correspond to the 
nonlinear and linear analysis of the discontinuous case for a — 0.7, respectively. 
The latter curve lies well below the former one, indicating that the linear analysis 
leads to an underestimation of the fluctuations. This is a general fact holding for 
all values of a in the meaningful range [0, 1]. The difference must be attributed 
to the sporadic amplifications due to the discontinuity: it is remarkable that 
the finite-size spectrum is independent of the initial condition. 

It is now important to test whether the difference between linear and non- 
linear curves persists also when the discontinuity is removed, i.e. in the CSM. 
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Figure 5: Probability distributions of the exponential growth rate (^2|) in the 
large time limit for the stochastoc model (10) in the 0-dimensional case, with 
a — 0.7. Lines 1 and 3 refer to the continuous case (A — 0.0097), with a finite 
size or an infinitesimal initial perturbation, respectively, while lines 2 and 4 refer 
to the discontinuous case (A — s- 0), still with a finite size or infinitesimal initial 
perturbation, respectively. Deviation from straight line of curve 2 is due to finite 
time effects. In the inset the difference between curves 1 and 3 is magnified. 



For the sake of simplicity, we fix A equal to a" and suppose that both n and 
no are non- negative integers, so that u{t) is defined on a discrete subset of the 
unit interval, {a"}„=o,i,... • This assumption does not affect the main conclu- 
sions while it allows one writing a simple recursive equation for the probability 
distribution Q{t,n): 

fi 

Q{t+ 1,0) = ^Q(i,A:)a'=+\ 

_ r (l-a")Q(t,n-l) + a'-'+iQ(t,n + n) < 7i < n 
" \ (l-a"+i)g(t,n-l) + a"+iQ(i,n + n) n>n ^ > 

Equations (pT|), with the initial conditions (5(0, n) = Sn.ng f^nd boundary con- 
ditions Q{t,n) = Q y t > n + no, can be numerically iterated to obtain the 
finite-size multifractal distribution in the continuous case. 

Since the development of analytical techniques to determine i?A is by far 
more complex than in the previous case, we have limited ourselves to determine 
the multifractal spectrum numerically. The linear and nonlinear spectra are 
reported in Fig. || for A = 0.0097 (curves 3 and 1, respectively). Let us first 
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notice that the finite-size spectrum is again independent of the initial condition 
u(0) and this makes Ha a well defined quantity: the efi'ect of ■u(O) is just to affect 
the convergence to the asymptotic spectrum. Furthermore, we can see that the 
nonlinear curve lies well above the standard multifractal spectrum, indicating 
that it is not the discontinuity to be responsible for the difference between the 
behaviour of infinitesimal and finite-size perturbations. Moreover, the overall 
closeness of curves 1 and 2 reveals that the removal of the discontinuity does 
not introduce significant differences in the finite-size spectrum. 

Finally, let us closely compare Ha (A) with /ia(A): the inset in Fig. || re- 
veals that the two coincide for A, A < Ac (see the dashed line), while above 
some critical value Ac the finite-size spectrum continues as a straight line, while 
Ha decreases faster. In the language of thermodynamic formalism, the linear 
behaviour of Ha is suggestive of a phase transition, from small values of A, 
that are correctly described by the linear analysis, to large values of A, where 
the finite character of the perturbation cannot be neglected. This point would 
certainly require a more detailed analysis to provide a more solid background to 
the above arguments, but we avoid this as it would drive us too far from what 
is the main goal of the present paper. 

The analysis performed in this section has allowed us introducing well defined 
observables to deal with finite-size perturbations. We cannot, however, avoid 
commenting on an alternative class of tools that have been devised to deal 
with this problem, although mainly in the context of truly chaotic systems. 
Since in realistic physical conditions perturbations are always finite, it is very 
tempting to introduce a growth rate to characterize perturbations of different 
sizes. However, serious conceptual problems are immediately encountered if 
one tries to define truly finite-size Lyapunov exponents. On the one hand, the 
finiteness of the size induces a dependence on the norm utilized, on the other 
hand it requires that the Lyapunov exponent has to be defined for a finite-time 
resolution, since perturbations change size over time. This latter implication 
is rather crucial in that so defined Lyapunov exponents are not self- averaging 
quantities (see Ref. jl6| for a detailed discussion of the problem). 

In spite of such limitations, finite-size Lyapunov exponents may carry useful 
information, although one has to be careful in interpreting them. Let us, for 
instance, look at 



where u' is the first iterate of u in the DSM and the average is performed over 
all possible realizations of the stochastic process. Simple algebra yields 



In the limit it — > 0, Ao(m) reduces correctly to the true Lyapunov exponent log a. 
If a is not too small (a > 1/e, where e is the Neper number), Aq(u) has a maxi- 
mum value for an intermediate value of u and, more interestingly, the maximum 




(22) 



Aq (u) = In a — aw In(au) 



(23) 
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value of A is larger than 0, if a > exp(— 1/e) = 0.692... Therefore, this result is 
suggestive of a phase transition from a regime (small a) where perturbations of 
all sizes decrease, to a regime where sufficiently large pertubations expand (and 
can, in principle, self-sustain). However, we already know that this conclusion 
is incorrect: the reason is precisely that Ao(u) is an average quantity, and fluc- 
tuations must be taken into account. In fact, a different scenario arises, if we 
define the contraction rate by taking the logarithm of the average expansion 
factor. In this case, we obtain 

Ai(m) = lna-Mn(2 - aw), (24) 

an expression that can be larger than zero even at u = 0, where it actually 
attains its maximum value! In practice, one has to be very careful in drawing 
meaningful conclusions from any expression of the finite-size Lyapunov expo- 
nent. In particular, since u can never be larger than 1, while it can become 
arbitrarily small, the typical negative contraction rate operating at small it- 
values eventually wins, making always u = the only stable fixed point. 



4 The spatially extended case 

The past study of stable chaos has revealed that no qualitative difference exists, 
as long as strictly finite chains are considered, with respect to standard stable 
systems. After an exponentially long transient, a periodic behaviour is always 
attained. As it was remarked in Section 11, a similar scenario is indeed exhibited 
by our stochastic models. 

Eq. (p^ ) shows that in the absence of any coupling, there is a finite prob- 
ability P(m(0)) for an initial perturbation m(0) to be contracted for an infinite 
amount of time steps and thus of being effectively absorbed. It is worth defining 
here 

UAiit) ~ maxwi(t) . (25) 

■i 

In particular, this allows us applying the above reasoning to the coupled case, 
with mm(0) playing the role of u{0); the probability P that an arbitrary, spa- 
tially extended, perturbation is contracted forever on every site i satisfies the 
inequality 



P{{U^}^ 



> 



Y[il-UMiO)a" 



.71=1 



p{uMm 



(26) 



where P{um{0))^ is small but finite for every finite lattice length L. Thus, 
the average time r needed for any perturbation to die out (i.e. to be contracted 
below some arbitrarily small thresold value) is finite and does not grow faster 
than exponentially with the system size. 



T < 



p{uMm 



(27) 
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Figure 6: The space-time and ensemble average value rh of the perturbation Ui{t) 
as a function of the contraction rate a. The dots refer to the result of numerical 
simulations with lattices of size L ~ 1000; the solid, dashed, and long-dashed 
curves refer to the mean field analysis, the Gaussian approximation discussed in 
the appendix B, and the numerical integration of the Frobenius-Perron equation 
(p9|), respectively. 



In the active, i.e. chaotic, phase where perturbations propagate with a finite 
velocity, r grows exactly exponentially (with the system size), since the only 
way for a perturbation to die out is to be contracted in all sites, while in the 
inactive phase, the latter is only a sufficent condition and not a necessary one, 
and perturbations die out on significantly shorter time scales. 

At variance with the 0-dimensional case, direct numerical simulations of the 
DSM in a 1-d lattice reveal the existence of a regime where finite perturbations 
self-sustain. This is clearly shown in Fig. ^, where we have reported the space- 
time and ensenble average value m of the perturbation for different values of the 
contraction rate. Above a = 0.60(5), nonzero amplitudes are actually observed. 
It is, therefore, crucial to understand the reason why the spatial interactions 
can stabilize finite perturbations in spite of the diffusive nature of the coupling. 

The natural extension of what we have learnt in zero-dimension consists 
in looking at the joint probability distribution R{t, ui, U2, . . . , Wn, . . .) over the 
whole lattice. If the system is sufficiently above the transition to stable chaos, it 
seems reasonable, in a first approximation, to neglect spatial correlations. This 
is certainly incorrect for those sites that are close to propagating fronts but the 
fraction of such lattice sites is definitely negligible. Accordingly, we approximate 
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the joint probability distribution as a product of single-site probabilities Q{t, u), 



R{t, ui, U2, ■ ■ • , Wn, ■ • •) ~ Q{t, ui)Q{t, U2) . .. Q{t, u„) (28) 

Within this approximation, the single-site probability distribution correspond- 
ing to the stochastic dynamics (|l^) satisfies the following Frobenius-Perron 
equation, 

Q{t + l,u) = g{u) Y[dv,Q{t,v,)d iu- —"^uA 

i=i \ i=i ) 

+ aA YidviQ{t,v^)S iu- ——^Vi\ 

.ex. .00 A' / a ^ \ 

Q{t + 1,1) = J duuj '[[dv,Q{t,v,)6 lu- —Y,^A ' (29) 



i=l 

where 6 is the Dirac's distribution, Uj is the amplitude of the perturbation in 
the «th neighbouring site, N is the number of democratically coupled sites (for 
later convenience we leave N unspecified - notice that TV = 3 in a 1-d lattice 
with nearest-neighbour coupling) and 



5(w) 



1-u A<u<l 
1 - aA < u < A 



It is easy to verify that the support of the single-site probability distribution 
Q{t,u) remains confined to the unit interval, provided that this holds true for 
the initial condition as well. Furthermore, due to the factorization hypothesis, 
the space and ensemble average m{t) defined in Section II coincides with the 
simplest ensemble average, i.e. the mean value of Q(t,u). 

In Fig. 1^ we have plotted the single site probability distributions obtained by 
directly iterating the stochastic model and the approximate Frobenius-Perron 
equation ( p9|) (solid and dashed line, respectively). The reasonable overlap 
confirms the validity of the factorization hypothesis (at least away from the 
critical region). The mean value, equal to 0.773 in former case, compares with 
0.770 in the latter one, while the variances are respectively equal to 0.063 and 
0.061. Such small differences are due to the different behavior of the probability 
distributions for small values of u. 

For iV = 1 (no coupling) , Eq. ( p9| ) corresponds to the 0-dimensional dynam- 
ics discussed in the previous section, and the evolution equation is exact. It 
reduces to the linear equation, 

Q{t + l,u) = g{u)-Q{t,-u) + aA'^Q{t,aAu) 
a a 

/•oo 

Q{t+ 1,1) = a / duuQ{t,u), (30) 
J A 
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Figure 7: Stationary, single-site probability distribution Q{u) for a — 0.7 and 
A = 0.01 computed by numerical simulation of the CSM (dashed line) and by 
numerical solution of Eq. ( [29| ) for = 3 (solid line) . Both distributions have 
been obtained subdividing the unit interval in 100 channels. The J-peak in 
u = I has been cut to magnify all other details. 



The only fixed point of this equation is Q{t, u) — (5(u), i.e. the absorbing state. 
The multifractal spectrum discussed in Sec. Ill is nothing but a sophisticate 
characterization of the convergence towards such a fixed point. 

On the opposite side of the 0-dimensional limit, there is the mean- field ap- 
proximation that corresponds to the limit N oo. In this limit, statistical 
fluctuations vanish and the dynamics reduces to the evolution of the mean value 
m{t), that reads 

I m{t)[2a^ahn(t)] if m(i) > A , 
+ " \ TO(t)[2a-a2A] ifm(t)<A. ^'^^> 

For a < ac = (1 — \/l — A)/A, Eq. (|l]) displays the stable fixed point mi = 0. 
This is the same regime found in 0-dimension and corresponds to the eventual 
absorption of any initial finite difference. Increasing a above Oc, the system 
undergoes a bifurcation: mi becomes unstable and a second (stable) fixed point 
m2 = (2a — l)/a^ appears. In the discontinuous limit (A — s- 0) Oc — > 1/2. In 
Fig. ^ the predictions of the mean-field approach are compared with the results 
of direct simulations: we see that, in spite of the approximations, there are 
no severe differences and the critical point is underestimated by approximately 
17%. 
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It is rather instructive to notice that the predictions of the mean-field anal- 
ysis do coincide with the finite-size Lyapunov exponent Ai (for the sake of 
simplicity, we limit ourselves to consider the DSM). Since the mean field ap- 
proximation reduces the CSM dynamics to the evolution of a single variable 
m{t), the comparison can be performed by interpreting ln[m{t + l)/m{t)] as Ai 
and recalling that Eq. (^J) has been derived for the DSM only. Even though 
Eq. ( p^ can be obtained by Eq. (|3^), there is an important difference between 
the consequences of the two results. As we discussed in Section III, fluctuations 
keep u* = stable for every value of a, while in the case of the mean field 
analysis, the lack of any fluctuation, due to the formally infinite number N of 
neighbours, implies that the fixed point m = is truly unstable when a > 1/2. 

The absence of fluctuations in the mean field limit implies that all defini- 
tions of the flnite-size Lyapunov exponent are equivalent. Therefore, we observe 
the same scenario previously observed for the standard maximum Lyapunov ex- 
ponent: the diffusive coupling shifts the Lyapunov exponent from the average 
value of the logarithm of the multiplier (the so-called quenched average holding 
for the single map) towards the logarithm of the average multiplier (annealed 
average, predicted by the mean field analysis). The important consequence of 
this shift is that, in the present context, it can change the stability of the u = 
solution leading to the onset of the chaotic phase. 

Anyway, one should not forget that the mean-field anlaysis provides an ap- 
proximate solution. For N finite and strictly larger than 1, Eq. ( p9| ) defines 
a non trivial evolution operator Q in a functional space. In practice, one can 
expand the evolution equation into an infinite set of equations for, e.g. the mo- 
menta Mfc of Q. An approximate solution can thus be found by either suitably 
truncating the hierarchy of equations or introducing a closure Ansatz. In App. ^ 
we parametrize the probability distribution as the sum of a (^-distribution and 
a Gaussian. This allows us deriving three evolution equations for the DSM. 

In both cases that we have investigated {N = 2, 3), we flnd a scenario similar 
to the one predicted by the mean field analysis. There exists a critical value Oc 
(equal to 0.548 for N — 2 and to 0.536 for iV = 3) below which the dynamics 
is characterized by the stable fixed point A — 1, v~Q, V — (corresponding 
to the absorbing state Q{t,v) = S{v)) and above which the previous solution 
becomes unstable, giving rise to a stable nontrivial solution. The dependence 
of m on a reported in Fig. ^ (see the dashed curve in Fig. ||) indicates that 
the critical value predicted by this analytic approach improves the mean-field 
estimate, but the growth of m above threshold is not as good as one would like. 
In fact, there is a qualitative difference with the mean- field approach: a further 
bifurcation (at a = a = 0.855 for iV = 2 and d = 0.869 for iV = 3), where the 
new solution destabilizes too. Such a bifurcation and the slow growth of m with 
a are both consequences of a defect of the approximation: the support of the 
Gaussian extends out of the unit interval. This unphysical property becomes 
increasingly important as soon as the average amplitude of u is comparable with 
1. 
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Indeed, a better agreement with the direct simulations is obtained by iterat- 
ing numerically the Frobenius- Perron equation (see Fig. |^). Neither simulations 
performed with A = nor with A = 0.01 reveal the second bifurcation found 
with the Gaussian approximation, confirming that it is an artifact of the ap- 
proximation. In the continuous case, the bifurcation occurs at the critical value 

= 0.591 ... (for = 2) and = 0.568 . . . (for N = 3), to be compared 
with the mean field prediction — 0.501 .... In the discontinuous case we find 
Oc = 0.585 . . . (for = 2) and = 0.567 . . . (for iV = 3), to be compared with 
the mean field prediction Oc = 0.5. 

From the data reported in Fig. |6| for = 3 we see that the factorization 
hyptohesis reproduces fairly well the behaviour of the full stochastic model ev- 
erywhere except for the transition region. This is not unexpected as it is well 
known that the correlation length diverges in the critical region. 

5 Open problems and conclusions 

In this paper we have shown that a simple stochastic model, specifically designed 
to simulate a different response to finite and infinitesimal perturbations, is able 
to capture the key features of irregular behaviour in linearly stable systems. 
In particular, we have seen that replacing the sequence of jumps generated 
by the CML dynamics with a genuine stochastic process allows for a faithful 
reconstruction of the front propagation. The main theoretical advantage of the 
stochastic model is the disentanglement between the generation of a pseudo- 
random pattern and the evolution of perturbations. In reality the two issues 
are interlaced: their separation has allowed us clarifying under which conditions 
(finite) perturbations can be effectively sustained throughout an infinite lattice. 
In particular, a full consistency exists between the CML and the stochastic 
model in the chaotic regime, since we can state that the amplification of finite 
perturbations contributes to sustain an irregular regime. 

On the other hand, the transition to the ordered phase observed in the 
CSM/DSM models does not reproduce the analogous behaviour displayed by the 
CML. In this latter case, it was observed that the critical region is not point-like, 
but rather extended to a what has been called "fuzzy region" , where ordered 
and chaotic dynamics alternate in a quite irregular manner . The reason for 
the difference is that in the CML model, the absence of local chaos makes the 
sequence of multipliers increasingly less random in the transition region. In the 
stochastic model, instead, the randomness of multipliers is always assumed a 
priori. In spite of such a difference, it is nevertheless instructive to notice that 
a transition persists without modifying the stochasticity in the real space. 

A more precise analogy for the transition investigated in the previous section 
is provided by the correspondence with the problem of synchronization in the 
presence of external noise. In fact, in this latter context, the noise represents 
the (unvaoidable) source of stachasticity in the synchronous as well as in the 
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asynchronous regime. Since a recent numerical study of the synchronization 
transition in hnearly stable system has suggested that it belongs to the univer- 
sality class of Directed Percolation (DP) [Q, it is tempting to verify whether 
the same holds true in our stochastic models. 

As already mentioned, even though the front velocity vp is a. good order 
parameter to characterize the transition, it is quite difhcult to obtain a reliable 
estimate of the critical value of the control parameter Uc from the vanishing of 
Vp. In fact, finite-size and transient effects combined with the existence of wild 
fluctuations prevent a careful analysis. A more efficient method amounts to 
measuring the dependence of the so-called absorption time r on the system size 
L. This is defined as the time required for the space averaged perturbation u{t) 
to become smaller than some very small, but finite threshold F. Fluctuations 
of r can be efficiently reduced by averarging over a sufficiently large ensemble 
of initial conditions. In the active phase, r is expected to diverge exponentially 
with L (stable chaotic regime), while in the absorbing phase, it should depend 
at most logarithmically on L. Only at the critical point, r exhibits a power law 
dependence 

r(L,ac)~L^ (32) 

where z = vijv^ is the so called dynamical exponent. We have performed 
numerical simulations for both the CSM (with A = 0.01) and DSM, averaging 
over 3000 realizations of the stochastic process and over randomly sampled 
initial conditions. In the CSM we find = 0.6055 . . ., with z — 1.58 ± 0.02; in 
the DSM we obtain Uc — 0.6065 . . . and z = 1.56 ± 0.06. The errors have been 
estimated as the maximum deviation from linearity in the log-log plot that has 
been used for extracting the scaling law (^2|) (see, for instance, Fig. ^, where r 
has been plotted versus L for A — 0.01 and different values of a). These results 
agree with the most accurate numerical estimates of the DP value, z — 1.5807 
II- 

We also measured the critical exponent 6 associated with the temporal decay 
of the density of active sites p{t), i.e. those sites where Ui{t) > F: at the critical 
point p{t) is characterized by the scaling law 

p{t,ac)^t-^, (33) 

where S = P/vi (as usual vi, and /3 are the critical exponents respectively as- 
sociated with space and time correlation lengths and with the order parameter). 
By averaging over 3000 realizations and choosing a sufficiently large value of L 
to get rid of finite-size corrections we have found 8 = 0.150 ± 0.01 for the CSM 
and 0.155 ± 0.005 for the DSM, to be compared with the DP value 5 = 0.1595 
HI- 

Altogether, our simulations support the hypothesis that the transition be- 
longs to the same universality class as DP. This may look as an almost trivial 
result, since local spreading is the only mechanism for the propagation of per- 
turbations (or, in different languages, active sites, infections). However, the 
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Figure 8: Log- log plot of the absorption time r as function of the system size 
L in the continuous case (A = 0.01) near the critical value for different control 
parameter values: a = 0.6051 (dashed line), a = 0.655 (solid line), and a = 
0.6061 (long-dashed line). 

whole problem is definitely more subtle, as an absorbing state cannot be iden- 
tified so clearly. In fact, we have already seen that the determination of the 
absorption time requires to fix a somehow arbitrary threshold F, and the same 
is true for the computation of the active sites at a given time. Even though 
we have found that our results are independent of the choice of F (provided 
that it is small enough), this may appear as a numerical trick: there is noth- 
ing like a true threshold, since no matter how small is a perturbation, there 
is always some finite probability that it gives rise to a burst: this is contrary 
to the existence of a truly absorbing state. On the other hand, the argument 
presented in the beginning of the previous section to convince the reader that in 
a finite chain any perturbation eventually dies out, confirms the existence of an 
absorbing state: since the smaller is um the more likely is that the perturbation 
keeps being absorbed, any perturbation ?i,:(t) has a finite probability to enter 
an infinite "contraction loop" in which every site i is monotonously contracted 
for any time larger than t. 

It would be nice to put our qualitative arguments on a more rigorous basis, 
by definining a suitable finite-size Lyapunov exponent that is negative below 
some threshold to indicate that u* = is a truly absorbing state. However, 
it is not clear whether this could be accomplished, since we know that in the 
chaotic phase, u* = should be at the same time "macroscopically" unstable 
since perturbations eventually drive the system towards the only stable state 
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and "microscopically" stable to mean that small enough perturbations have to 
be absorbed. In the future, we hope to be able to clarify whether it is possible 
to define an indicator that contains both messages. 
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A Multifractal distribution 

In this appendix we report the analytical calculation of the (constant) slope of 
the multifractal distribution (|l^) for the zero dimensional DSM. We shall prove 
that H'{A, u(0)) (the prime denotes derivative with respect to A) is independent 
of A and increases monotonously with the contraction rate a from -1 (for a — *■ 0) 
to (for a 1). 

First of all, note that the condition < n < uq + t, together with Eq. (|l8|), 
implies that the support of H{A) is confined in the interval [lna,0], since the 
ratio p = n/t can take values between and 1 and in the t —t oo limit no/t 
vanishes. 

We are now interested in the A > Ina case (i.e. the one in which the system 
received at least one kick), where, as it was stated in Section III, the actual size 
of the perturbation u at time t is unambiguously determined by the time n <t 
elapsed since the last kick (see Eq. (p^) and thus, for any finite i, both u and 
A can assume only a discrete set of values labelled by n. 

Let us now denote the discrete A-derivative of a generic function /(A) as 

Z,.,/,./<A±^M, (34) 

From Eq. ( |l3|) , we can approximate the derivative of the multifractal distribution 
H as 

i7'(A,.(0))=limM!iZ(MM!^, (35) 

where AA is naturally fixed by Eq. ( p^ ) and the discrete character of n = 
Inw/ Ina, 

AA = A,+i (t, 7i(0)) - A, {t, ,i(0)) = . (36) 
As, for t — > oo, AA goes to 0, Eq. (pq) becomes asymptotically exact. 
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Morover, Eq. 
ities writing 



allows one switching to the Q rapresentation of probabil- 



1 



iJ'(A,u(0)) = lim In 
In a t 



■P(A + AA,i,M(0)) 



In a t- 



lim In 



Q(i,n + 1) 



Q{t,n) 



(37) 

(remember also that u{0) = a"°). Making use of Eqs. ( p^ and (pi]), we obtain 



H'{A,u{0)) 



1 

In a 



In 



/ G(t-n~ l,no) 
V G{t-n,no) 



In (1 - a 



ri+l 



(38) 



In the limit t — > oo, the last term in the r.h.s. vanishes, provided n diverges 
too (i.e. A < 0). We are thus left with the following equality, which holds true 
for any A in the open interval (In, 0), 



where 



H'(A,um = - 
r]{a, no) = lim 



In 77(0, no) 
In a ' 

G(t + l,no) 
G(t,no) ' 



(39) 



(40) 



the n-dependence has been eliminated by shifting the time origin and the de- 
pendence on the parameter a has been made explicit. From this equation, we 
see that the problem of determining the slope of the multifractal spectrum is 



equivalent to an eigenvalue problem. In fact, we can formally write Eq. (20) as 

v{t,no)^Mv{t-l,72o), (41) 

where v{t,n()) is the infinite-dimensional vector 

v{t, no) = {G{t, no), G(t - 1, no), G(t - 2, no), . . . , ) (42) 

and M is a linear infinite-dimensional operator. Therefore, 77(a) is nothing but 
the maximum eigenvalue of the operator M. 

From the explicit expression of the recursive equation for G{t, no) (see Eq. ( |2C| ) 

t-i 

G{t, no) = G{t-1, no)P(0, 1) a+G{t-2, no)P(l, 1) a^+. . .-ha"° J| (l-a'=+"°)a*, 

k=l 

(43) 

we see that, since both G{t,no) and P{t,l) are not larger than 1 (they are 
probabilities), the sum in the r.h.s. is bounded from above by the sum of 
the first t powers of a. In the limit t ~f 00, the function G'(i, no) converges 
exponentially fast, so that we are able to truncate Eq. (^) to order k with an 
arbitrary precision. This truncation makes the problem numerically solvable. 
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as the operator M can be approximated by the finite dimensional matrix Mk 
(here and below, the subscript k stands for fcth order approximation) that we 
report here below, 



Mi. 





a2 


as ■ 


• Oik 


1 





• 


■ 





1 


• 


• 








1 • 


• 








• 


• 1 



(44) 



1 / 



where aj = P{j — 1, 1). Numerical estimates of rik{a) indicate that fc = 10 
suffices to attain a good convergence in the whole range of a values between 
and 1. For k = 2 the eigenvalue ri2{a) can be computed analytically yielding. 



g^(A,^(0))^-l-^ln ^ + f^ , 
ma 2 



(45) 



which goes monotonously from -1 (a ^ 0) to (a ^ 1). 



B An approximate solution of the 1-dimensional 
problem 

The Frobenius-Perron operator defined by Eq. ( p9| ) is a functional equation 
whose solution requires projecting it onto a finite dimensional space either via 
some truncation or a suitable closure hypothesis. Looking at a typical shape 
of the stationary probability distribution in Fig. ^ we can see that the multi- 
peaked structure slows down the convergence of an expansion in either moments 
or cumulants. Accordingly, we have preferred to approximate the probability 
distribution as the sum of a Dirac's S distribution centered in 1 (the reinjection 
point) plus a Gaussian distribution, centered around a point to be determined 
self-consistently, 



Q{t, v) = AtSiv - 1) + i-=^ cxp 



2Vt 



(46) 



We can see that Q{t,v) is parametrized by three quantities: At, the probability 
of the 6 component, Vt the variance of the Gaussian, and Vt, its average value. 
The mean value m{t) of Q{t,u) is therefore equal to 

m{t)^l-At + AtVf (47) 

It should be noted that our approximation is formally "unphysical" since the 
support of any Gaussian function is not restricted to the unit interval, but we 
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expect our Ansatz to be reasonably correct as long as the probability to be out 
of the unit interval is small enough. 

Entering the above definition of Q{t,v) into Eq. ( [29| ) and computing sepa- 
rately the new weight of the 5 component, the average and the variance of the 
Gaussian, we obtain three evolution equations. In the simplest nontrivial case, 
N — 2, we have 



At+i 

Vt+l 

Vt+i 



a[At{l-vt)+vt] 

-v^ [2 - At] y + Vta [1 - aAt] - a 



a , . , 1 — a 

l + -[V-2 + At)-^^ 



-At\a - -a 



Vt 



^^-{l-af 
3 



a + Vt 



At{2-a)-a^ 



a^ — a 



1 



-At{l - Vt)-y - -a-'Vt + a\Vt - 2) + 2a - 



2a(l - a) 
I -At 



At[a{Vt-l) + 2]^ + a^ + ^{Vt-2) 



I -At 



For = 3, the case that we have investigated in detail as it corresponds to a 
1-d lattice of democratically coupled maps, the equations read 



At+i 

Vt+l 

Vt+i 



a[At(l ~ vt)+vt] 

a^ 4 
-v'i [3 ~At]—+ vt{a - 3 a^^t) - a 



l + -iVt-3-2At) 



1-a 



I -At 



-^A^to" + ^^t(5a - 2)a - (1 - af 



-A 



2„3 



\At{2a - 3)a^ + a^ - a + 



2vt_ 
3 



2„3 



Afa- 



At{-?, + 2a + Vta) - 2a -a^ - 1 + —(2 - 3a)a 



1-a 
l^t 

2 



2 



1 - 



A? f ^a^ + a^ - a 1 - At 



—a 
9 



Vt—-a 



2{a-a^) 
l-At ' 



In spite of the great simplifications involved in the derivation of both sets of 
equations, it is still impossible to obtain an analytic expression even for the 
critical point Oc. The results of the numerical solution of Eq. ( p9| ) and (^) are 
discussed in the text (see Sec. IV). 
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